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We address a fundamental issue in quantum mechanics and quantuin information theory, the 
generation of an entangled pair of qubits that interact solely through a third, semiclassical degree 
of freedom, in the framework of cavity quantum electrodynamics. We show that finite, though 
not maximal, entanglement is obtainable in the classical limit, at the price of a diverging effective 
interaction time. The optimal atomic entanglement derives from a trade-off between the atomic 
entanglement in a sub- wave packet and the purity of the atomic state. Decoherence by photon loss 
sets an upper limit on the degree of excitation of the cavity mode, beyond which the achievable 
entanglement decreases as the inverse mean photon number to the sixth power. 
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I. INTRODUCTION 

Field-atom entanglement is one of the hallmarks of 
strongly interaeting cavity quantum electrodynamics 
(CQED) systems. This fundamental process is often used 
as a building block in the formation of entangled states 
of two or more atoms that serve as a resource for quan- 
tum information processing. A notable example is the 
generation of a two-atom Bell state by means of consec- 
utive interaction with a resonant microwave cavity mode 
prepared in the vacuum state |l|. In effect, the cavity 
mode stores the entanglement with one atom, that is 
then transferred to the second atom. The question ad- 
dressed in this work is: Can entanglement be generated 
when the initial field state becomes semiclassical? On one 
hand quantum fluctuations persist also in the semiclas- 
sical limit, while on the other hand the correspondence 
principle states that in this limit classical physics should 
be recovered, where the notion of entanglement is not 
even defined. 

The semiclassical states of the field mode are wave 
packets that are localized in phase space, where the 
canonical coordinates are the field quadratures, such as 
coherent states and squeezed coherent states, and have 
a large mean photon number. When a two-level atom 
interacts with a coherent state wave packet, its polar- 
ization undergoes Rabi oscillations whose amplitude ex- 
hibits collapse-revival dynamics [2[ as a consequence of 
the splitting of the initial wave packet in two mutually 
orthogonal sub-wave packets, and of their periodic col- 
lisions [3, |j|- Each sub wave-packet is a product of 
squeezed field state and an atomic state. After split- 
ting, the atom-field system is highly entangled, but slow 
atomic state evolution turns the field-atom state again 
into a product state 0, |j] , and the field state is a super- 
position of two well-separated wave packets, sometimes 
called a Schrodinger cat state. This semiclassical dynam- 
ics can be described as a fiow in a double phase space Q ■ 

Here we study a field mode interacting consecutively 
with two two-level atoms. The interaction of the field 
wave packet with the first atom causes it to split, as ex- 
plained, and to become entangled with the atom. When 
the interaction with the second atom begins, each of these 



sub-wave packets splits again, and the atom-field-atom 
system becomes a superposition of four wave packets (see 
Sec. ini Fig. [T]). Atom-atom entanglement can only be 
obtained if two or more of the sub-wavepackets overlap, 
and this occurs when the normalized interaction times of 
the field with the two atoms are equal. The system state 
is then in general a superposition of three wave packets, 
and the atom-atom state is an entangled mixed state. 

We show in Sec. IIIII that it is possible to generate en- 
tanglement that reaches a finite, less than maximal, limit 
when the photon number tends to infinity, although the 
time required to generate the entangled state diverges in 
this limit. The optimal interaction time is determined as 
a trade-off between the mixing of three pure atomic states 
and the degree of entanglement of one of them, as shown 
in Fig.[3l We consider the effects of photon loss in Scc. lIVI 
and show that the entanglement generation is degraded 
by dephasing of the sub wave packets, and that a coher- 
ent superposition of wave packets cannot be maintained 
when the mean photon number is larger than the inverse 
decoherence rate to the | power. Beyond this point en- 
tanglement generation proceeds only through small quan- 
tum fluctuation and it decreases as the inverse of the 
mean photon number to the sixth power, see Fig. |4l 

The analysis demonstrates that a semiclassical prepa- 
ration can serve as mediator for appreciable entangle- 
ment of atoms, if the decoherence is weak enough. For 
this purpose it is necessary that the field becomes entan- 
gled with the atoms and evolve into macroscopic superpo- 
sition states. In contrast, when decoherence is too strong 
to allow wave packet splitting, atom-atom entanglement 
is generated only through small quantum fiuctuations, 
and it is therefore weak. 

The resonant interaction of a field mode in excited 
states with two atoms has been investigated theoreti- 
cally in p-Q for simultaneous interaction and [3, llfl for 
consecutive interaction, revealing wavepacket splitting, 
slow atomic state evolution, and collapse-revival dynam- 
ics of entanglement similar to those observed in one atom- 
field system. The off-resonant interaction of a field wave 
packet with two atoms was studied experimentally and 
theoretically in [Tl], [l2l as a means of measurement of 
the cavity decay rate and its effect on the dephasing of 
Schrodinger cat states. The steady state entanglement 



of two atoms in a highly excited leaky cavity was stud- 
ied [13]. Our focus is the process of entanglement of 
two two-level systems by interaction with an intermedi- 
ary that tends to the semiclassical limit, for which CQED 
provides a prominent realization. 

Field-mediated atom-atom entanglement is often car- 
ried out in the highly-detuned Raman regime [ij, |l5| 
where the field is only virtually excited to avoid field- 
induced decoherence. In a sense this is the opposite limit 
to the one studied here, and indeed we show that the 
atom-atom entanglement is very sensitive to dephasing 
of the field wave packet. Furthermore, although it is of- 
ten argued that in the Raman regime entanglement is 
independent of the field state, the Raman regime entails 
an arbitrarily large detuning when the field state photon 
number tends to infinity. 



II. ATOM-FIELD-ATOM INTERACTION 

We study the dynamics of two two-level atoms inter- 
acting consecutively on resonance with a single electro- 
magnetic mode. Modelling the interaction by the Jaynes- 
Cummings Hamiltonian in the rotating wave approxima- 
tion, the Hamiltonian is H{t) = Hi, < t < ti, and 
H2, ti < t < ti + 12, where 



Hk 



fujj{a' a + (7\(7k) + hVtk{a) (7k + a(7\) (1) 



where to is the frequency of the electromagnetic mode 
with energy states |n), n = 0, 1, . . . created by a\ and 
coherent states |a), ftw is the level spacing of the two- 
level atoms with energy states \g)^. , |e)j, (atom k) raised 
by aj,, and hflk are the field-atom interaction energies. 
The system is prepared in a product state |a)®|(7)^®|g)2. 
When the field interacts with the first atom, the en- 
ergy states are polariton states |n)_|_ = -y=[\n— l)(8)|e)^± 

|n)(g) |g)^), with energies H{nijj ± ^ynQl) , and the absolute 
ground state \0) (^\g)-^. The field-atom state evolves into 
a superposition of two products of a squeezed coherent 
state and an atomic state [2|-lj]. Expressing the Hamil- 
tonian in terms of polariton number operator h, and the 
projections P± on the ± subspaces 13,^^1= H+P+ + 
H-P-, where H± = fi{nuj ± ynf^i), the initial state is 
■^i\a)^ - \a)_), where |a)± ^ ±%/2P±(|a) ® \g)^), and 
the field-first atom state at time ii is 
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For large |a| the dynamics generated by H± is well- 
approximated by classical dynamics in the polariton 
phase spaces, that is, phase spaces corresponding to the 
± sub-Hilbert spaces, where the canonical coordinates 
q,p are the field quadratures. For convenience, we join a 
|0)_|_ state to the ± subspaces (respectively) with occupa- 
tions that remain exponentially small throughout. The 
initial states in both phase spaces are then coherent state 
Gaussian wave packets, and they evolve according to the 



rules of semiclassical phase space dynamics [l6[ . In this 
approximation the wave packet evolution is determined 
completely by classical data generated by the classical 
Hamiltonians 



H. 
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^{q'^ + p' + h)LO + ^^{q^+p^ + h)ni (3) 



obtained from the Weyl phase-space representation of the 

(cl) 

quantum Hamiltonians. H)^ generate nonlinear oscilla- 
tions, that is rotation in phase space with an amplitude- 
dependent frequency. 

The wave packet propagation consists of three parts: A 
phase-space translation from the initial point {q,p), a = 
-^{q + ip) to the final point (g±,p±), a± = -^{q±+ip±) 
determined by the classical trajectory of the center of the 
wave packet, squeezing determined by the phase-space 
deformation generated by the nonlinear oscillations with 
squeeze parameters ^±, and an overall phase factor e~"^± 
determined by the classical action of the classical orbits, 
so that the wave packets evolve to the squeezed coherent 
state 
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b being the polariton annihilation operator defined with 
the usual properties b\n)^ = -yn |n — 1)_|_, [b,b^] = 1. 
The values of the wave packet parameters in a frame 
rotating with angular speed uj are (f>± = ^(±71), a± = 
"(±7i): C± = ^(±71) with 

0(7) = i|ap-iarctan(i7) (6) 

a(7) = ae~3'' (7) 

^(7) = arcsinh(i7)e*('^+^'''='=°*(37)) (g) 

Here 71 = ^^ is twice the phase space angle of rota- 
tion, clockwise and counter-clockwise (respectively), of 
the + and — wave packet; ^i^ is the classical (angular) 
frequency of nonlinear oscillations. 

The atom-field states can be expressed in terms of 
photon states with the help of the identity \n)^ — 
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where 70 = arga. Here we used the fact that the 
squeezed state wavepackets are localized in phase-space, 
so that they are approximate eigenstates of a and a^ with 
eigenvalues a(±7i) and a(±7i)* (respectively). 

The atom- field state |I) is therefore a superposition of 
two wavepackets that are well-separated in phase space 
except when 71 is close to an integer multiple of 27r. The 





FIG. 1: The Wigner function of the field state Tratoma |2) (2|, for unequal effective interaction times 71 7^ 72 (left) and equal 
interaction times (right), numerically calculated for a = 4. In the left panel the arrows signify the direction of phase space 
rotation of the four sub wave packets. The Wigner function of the initial coherent state wave packet is circularly symmetric, 
and the deformation of the four sub wave packets in the left panel, and the two flank sub-wave packets in the right panel is a 
consequence of the squeezing. 



splitting of the wave packet generates field-atom entan- 
glement that is almost maximal when 71 > |q^|^^ and 
decreases to zero for 71 = tt Q. 

The interaction with the second atom proceeds anal- 
ogously. The final field-two atom state at time ti + ^2 
is 
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where each of the two wave packets |1)_|_ 1.9)2 splits 
again to two sub wavepackets, that continue to undergo 
phase-space rotation and squeezing according to the field- 
second atom polariton sign. |2) is therefore a superpo- 
sition of four localized wave packets labeled by two sign 
choices 
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where F^^ = {r^n + S72) and 72 = ^. 
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III. ATOM-ATOM ENTANGLEMENT 

The field- two atom state |2) is a superposition of four 
sub- wave packets each with a well-defined atomic product 
state. Superpositions of product states can give rise to 
entanglement, but for most values of 71 and 72 the four 
sub-wave packets are separate in phase space (see Fig. 
[T]), and label the different atomic components, so that 



the two- atom density matrix pa = Trfidd |2) (2| is a clas- 
sical mixture of four product states in the limit |a| ^ 1, 
which is by definition separable. The atomic states form 
a coherent superposition when the corresponding sub- 
wave packets overlap in phase space. We therefore let 
7 so that the field factors in the |2_| ) and 



71 = 72 

1 2 (-) terms are identical and equal to the initial coher- 

(10) ent state |a), while being separate from the field factors 

of |2++) and |2 ), unless 7 is an integer multiple of 27r. 

If 7 is not an integer multiple of tt, the field factors of 

|2++) and |2 ) are also phase-space separate, so that 

the atomic state is a mixture 



p, = i(|0(Z|+2|c)(c| + |r)(r|) 



(13) 



of the states 
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+ ie'T" sin(i7) |e)i ® \g)^ -\9)i® [5)2) ^ and 
V) -^(Iff)i + e'^oe-t'' |e)i) ® {\g)^ + e^-^^e"^^ le)^) . 

The state |c) is entangled for all 7 not an integer mul- 
tiple of TT, with monotonically decreasing entanglement 
entropy as a function of 7 between and tt in the limit 
I a I —J- 00 (see Fig. [2]), while the states \l) and \r) are sep- 
arable. The full two-atom state is entangled if pa can- 
not be expressed as a mixture of separable states. The 
entanglement of formation jl7| is an entanglement mea- 
sure that vanishes if and only if the state is separable, 
displayed in Fig. [3] that shows that the two atoms are 
entangled for most values of 7 for finite a, and for all 7 
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FIG. 2: The entanglement entropy Ec of the center wavepacket two-atom state \c) (left), and the purity P of the full two-atom 
state pa (right) as a function of the phase-space rotation angle 7, shown for a = 4 (thin red), a = 8 (medium green), and 
asymptotically for |Qf| — > 00 (thick blue). The initial rise and later drop in the degree of entanglement of the full two-atom state 
can be understood as a trade-off between the increase in P and the decrease of Ec- P{a — >■ 00) is discontinuous at 7 = 0,7r, 
and these points are excluded from the graph. 
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FIG. 3: Left: The entanglement of formation Ef of the two-atom state pa as a function of the phase-space rotation angle 7, 
shown for a — i (thin red), a = 8 (medium green), and asymptotically for \a\ — s- 00 (thick blue). Right: The entanglement of 
formation (upper, blue) and negativity (lower, violet) as a function of 7 for \a\ — ^ 00. The two entanglement monotones attain 
their maxima (shown by vertical dashed lines) at close but different values of 7, showing that they are nonequivalent. 



except integer multiples of tt in the limit |a| — > cx3. When 
nonzero, the entanglement of formation of the two-atom 
state, being equal to the minimal weighted average en- 
tanglement entropy of pure states mixed in pa, obtains 
its maximum of « 0.21 at a kink singularity at 7 = ^. 
Unfortunately, there exist nonequivalent entanglement 
monotones for mixed states [l8|. As an example, the 
negativity, that is minus the sum of the negative eigenval- 
ues of the partial transpose of pa , obtains its maximum 
at 7 « 0.51l7r (see Fig. [3]). This result indicates that 
although there is no precise meaning for optimal interac- 
tion time, values of 7 near ^n are best for entanglement 
generation in this system. Furthermore, changing the 
initial condition of the second atom to | e) leaves the en- 
tanglement of formation unchanged, and slightly shifts 
the maximum of the negativity to 7 « 0.5137r. 

An important aspect of the two-atom entanglement 
generation process is that although the atom-field en- 
tanglement, and therefore the entanglement entropy of 



the center wave packet | c) , reaches maximal value after a 
short effective interaction time 7 of 0(|a|~^), the degree 
of entanglement of the full two-atom state is very low for 
such small 7. This is a result of the incoherent mixing of 
\c) with the flank wave packets |^) , \r) whose overlap with 
\c) is low for small 7, as follows from the low purity of pa 
for these 7, see Fig. [51 When 7 increases, the entangle- 
ment entropy of \c) decreases until it becomes separable 
for 7 = TT, but at the same time the overlap of |c) with 
1^) and I?'), and consequently the purity of pa, increase. 
As a result, the degree of entanglement of pa increases 
initially, reaches a maximum near 7 = ^tt, and decreases 
to for 7 = TT, as shown in Fig. [31 Although Eq. p^ 
is not valid for 7 an integer multiple of tt, the preceding 
argument is valid for these values of 7 showing that the 
entanglement of formation is a continuous function of 7. 



IV. EFFECTS OF DECOHERENCE 

The preceding results apparently contradict the corre- 
spondence principle by displaying entanglement, a purely 
quantum phenomenon, at the classical limit. However, 
a given value of entanglement is reached for a fixed 
^ " IcH^ ~ 1^' ^"^^ therefore the required interac- 
tion time diverges when \a\ — ^ oo. This is a demonstra- 
tion of the singular nature of the semiclassical limit: for 
a fixed interaction time classical physics is recovered as 
\a\ — >■ cxD, but the classical and the long-time limits do not 
commute. A similar phenomenon was observed in [19| . 
where the classical limit and the limit of large squeezing 
do not commute. 

Physically, longer interaction times imply a stronger 
effect of environment coupling, and it is this effect that 
guarantees the emergence of classical physics for large |a|, 
both because of stronger dephasing and because of longer 
interaction times. We demonstrate this statement using 
the standard Markovian model of cavity loss, with Lind- 
blad generator a and rate A, so that the master equation 
for the full system density matrix p is [20[ 

dtp^ -Jl[H,p\ + \\{a^ ap - 2apa^ + paU) (14) 

A superposition of distinct coherent states c(|a) -I- |/3)) 
experiences as a result of cavity loss, in addition to the 
overall decay with rate A, a much faster dephasing with 
rate Aja — /3p that affects the coefficients of the coherence 
terms in the density matrix proportional to \a) (/3| and 
its conjugate j2g |. 

In order to analyze this process in conjunction with 
the wave packet dynamics we assume that \\a\ <^ Sli, 1^2 
so that we can ignore the energy decay of the state. The 
full-system density matrix after the interaction with the 
first atom is then 

Pl = |1 + )(1+| + Xl |1 + )(1_| +.T* |1_)(1+| + |1_)(1_| 

(15) 
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with Xk = e '^fc ^' ^ ".A similar dephasing af- 

fects the coefficient of the twelve coherence terms during 
the interaction with the second atom. 

The significance of dephasing for the generation of two- 
atom entanglement is the reduction of the coherence be- 
tween the two wavepackets |2_| ) and |2 |_) that gener- 
ates the entangled atomic state |c) when they collide, so 
that the term |2^ ) (2 |_| and its conjugate in the full- 
state density matrix arc multiplied by xia;2 and (CC1X2)* 
(respectively). It follows that the term ^ \c) (c| in pa is 
replaced by 



\{\ca) {ca\ + a;ia;2 \ca) {cb\ + (a;ia;2)* \cb) {ca\ 



Pc- i 

+ \cb) (cb 



(16) 



with 



where upper (lower) signs correspond to A (B), respec- 
tively, while the terms proportional to |/) {l\ and \r) {r\ 
are negligibly affected by decoherence. For |xia;2| < 1 Pc 
is the density matrix of a mixed state. Since \ca) and \cb) 
are product states, the degree of entanglement of pc, and 
therefore p^ , is a decreasing function of the decoherence 
rate A. 

It follows from Eq. (fT6|) that the degree of entangle- 
ment depend on the decoherence rate only through the 



combination y 



2A|ct|-' 



Fig. |4] shows the entanglement 
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of formation as a function of y and 7. Evidently, the max- 
imum achievable entanglement decreases with increasing 
y, and this maximum is achieved earlier. Detailed analy- 
sis shows that when y is large the maximum of the entan- 
glement is achieved at 7 = -, and its value is proportional 
to y^logy. This phenomenon has the simple interpreta- 
tion that a shorter effective interaction time allows less 
time for entanglement generation, but also less time for 
decoherence, and that this trade-off leads to a shorter 
optimal interaction time for stronger decoherence; for 
large y the entanglement is completely destroyed when 
7 = J-. In this limit the interaction stops before the 
initial wave packet splits so that entanglement is gener- 
ated only by weak quantum fluctuations, and this allows 
a power law rather than an exponential decay in the de- 
gree of entanglement — the entanglement of formation for 
the optimal interaction time is approximately ^ of the 
entanglement achieved in an ideal system for the same 
interaction time, sec Fig. [5] (right panel). 



V. CONCLUSIONS 

Our first conclusion is that two qubits can be entan- 
gled solely by interaction mediated by a macroscopic 
field preparation, and the degree of entanglement tends 
to a positive value in the classical limit. The entan- 
glement process proceeds through splitting of the field 
wave packet into four components each carrying a differ- 
ent atomic product state and a subsequent merging of 
two of the sub wave packets to form a entangled super- 
position of two atomic product states. Thus, although 
the initial field state is semiclassical with a large pho- 
ton number and well-defined quadratures, the field nec- 
essarily evolves into highly non-classical states during the 
entanglement generation process. Nevertheless, semiclas- 
sical wave packet dynamics approximation for the prop- 
agation of the wave packet stays valid during the full 
entanglement generation process, although in a multiple 
phase space — one phase space copy for each sub wave 
packet. 

The splitting of the wave packet is a slow process where 
a macroscopic field state evolves by interaction with two 
qubits. Naturally, it is not an efficient method of creat- 
ing entangled pairs, as the interaction time required to 
obtain entanglement diverges in the classical limit. The 
correspondence principle is therefore nonetheless obeyed 
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FIG. 4: The entanglement of formation Ef of the two-atom state as a function of the phase-space angle 7 for y = ^ ^^ = 
0,0.15,0.4,0.7, 1 respectively (left) and y = 0,4,5,6,7,8 (right) from top to bottom, where A is the cavity loss rate and ili,2 
are the atom-field interactions. The dashed line in the right panel traces the majdma of Ef over 7 for given values of y. 



in the sense that classical physics, i.e. no entanglement, 
is obtained in the limit of large photon number for a fixed 
interaction time. 

The final two-atom state is a coherent superposition 
of two product states associated with the center wave 
packet, incoherently mixed with two additional products 
associated with the flank wave packets. Thus, the two- 
atom entanglement is determined by the overlap of the 
evolving atomic states. In particular, wave packet split- 
ting and atom-field entanglement are necessary but not 
sufficient for the generation of atom-atom entanglement. 

A subtler issue is the fact that the two atoms are not 
maximally entangled by interaction with the semiclassi- 



cal field mode. Although for short interaction periods the 
center wave packet has an almost Bell atomic state, the 
atom-atom entanglement is degraded almost completely 
by mixture with the flank wave packets, while for long in- 
teraction times the field and atoms decouple, the atomic 
state at this stage is separable. Unlike the vacuum field, 
therefore, the semiclassical field does not function as a 
quantum gate, and we conjecture that this is an example 
of a general principle. 
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